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Generating maximal entanglement between non-interacting atoms 
by collective decay and symmetry breaking 
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A simple scheme is presented for achieving effectively maximal pure-state entanglement between 
non-interacting atoms through purely collective decay and controlled symmetry breaking. The 
scheme requires no measurements or feedback or even knowledge of the initial states of the atoms. 
It relies on breaking the symmetry of the system Hamiltonian to ensure the existence of a unique 
attractive steady state and minimal control to achieve almost perfect overlap of this steady state 
with the maximally entangled singlet state. We demonstrate how our scheme can be implemented for 
two qubits encoded in hyperfine levels of atoms such as Rubidium in a lossy microwave cavity using 
only small magnetic field gradient. Error analysis suggests considerable robustness with regard to 
many imperfections including atomic decay, asymmetric atom-cavity coupling and frequency offsets. 
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Entanglement is an important resource in quantum 
physics, especially quantum information, and many 
schemes to generate entanglement have been proposed, 
using various tools including coherent control, mea- 
surements and feedback. In coherent control schemes 
(e.g., [1 ) dissipation tends to destroy entanglement cre- 
ated by Hamiltonian evolution and reduce the control 
effectiveness, but it can also be a resource for creating 
entanglement. Some interacting systems have entangled 
ground states and relaxation processes can drive the sys- 
tem into such states, known as cooling into entangled 
states [2 . More generally, careful design of the Hamilto- 
nian and interaction with the environment, i.e., Hamilto- 
nian and reservoir engineering, can ensure decay to a cer- 
tain state [SHS]. Reservoir engineering by measurements 
has been used to drive systems into entangled states [6 , 
and schemes relying on direct or indirect feedback have 
been proposed [7 . Entanglement induced purely by in- 
teraction with an environment has also been explored [8 , 
but the entanglement achieved this way is usually low- 
grade mixed-state entanglement and often transient, i.e., 
the system relaxes to a separable steady state eventually. 

Ideally, we would like to engineer systems that relax to 
a stable, maximally entangled, pure steady state regard- 
less of the initial state. For open systems described by 
Lindblad master equation this can be achieved by engi- 
neering the dynamics to render the steady state unique 
as uniqueness implies that any initial state "collapses" to 
this state [HIS]. Based on this idea, we show that effec- 
tively maximal pure-state entanglement can be created 
for atoms in a damped cavity driven by a classical field, 
even if there is no interaction between the atoms and no 
measurements are performed, i.e., purely by dissipation^ 
and the initial state is unknown. Of course, the singlet 
state for two identical atoms in a cavity is maximally 
entangled, but it is not attractive, i.e., if the system is 
initially in any other state, it will not evolve into the 



Gradient magnetic field 



Driving field 




Heavily damped 



FIG. 1: Two trapped atoms in a heavily damped cavity. The 
cavity field is slaved to the driving field. The quantization 
axis and field polarization direction are defined by a static 
magnetic field along z-axis. A small magnetic field gradient 
along X-axis permits variation of the atomic energy levels. 



latter. However, breaking the symmetry in the model, 
e.g., by varying the energy level splitting of the atoms 
slightly, results in a unique attractive steady state, and 
surprisingly, this is a pure state whose concurrence can 
be made arbitrarily close to 1 by suitable choice of the 
parameters, and which can be reached in a short time 
from an unknown initial state simply by adiabatically 
decreasing the induced asymmetry. 

There are many possible physical realizations for the 
type of system described but to fix the ideas and demon- 
strate the feasibility of the scheme we consider two 
trapped ultracold atoms such as ^^Rb, situated in a 
damped microwave cavity driven by a classical field as 
shown in Fig. [T] as a particularly promising implemen- 
tation. A magnetic field is applied to provide a quan- 
tization axis and photon polarization direction, and the 
qubit states are encoded in two hyperfine magnetic sub- 
levels such as It) = \F = l,mi? = 1) and ||) = \F = 
2^mF = 2). These basis states are stable with very long 
coherence times, i.e., once the atoms are prepared in 
a state such as the singlet state, we can expect them 
to remain in this state for a long time. The system 
could be prepared by loading the ultracold atoms into 



optical lattices to ensure the atoms are well localized 
and stationary. A detuning of the atomic transition fre- 
quencies of ±Auj from the cavity field frequency u can 
be achieved via the Zeeman effect by applying a small 
magnetic field gradient. The total Hamiltonian of the 
system is Htot = Ha ^ He ^ Haa + Hac + Hdc, where 
the Hamiltonians for the two atoms and the cavity are 
Ha = h{ujialai -\- uj2cr2cr2) ^ He = hujb^b and an (n = 1, 2) 
and b are the annihilation operators for the nth atom and 
intra-cavity field, respectively. In a typical setup with 
neutral atoms at least several microns apart direct inter- 
actions are negligible, Haa = 0- The atom-cavity interac- 
tion is given by Hac = h[gi{ai-\-al)-\-g2{cr2-\-crl)]0{b+b^) 
and the classical coherent driving field f{t) = Acos{ujt) 
adds the term H^c = iAdcos{ujt){b — b^)^ where d is the 
coupling strength. Including spontaneous emission and 
decay of the cavity field, the evolution of the total density 
operator pt of the driven atom-cavity system is given by 

Pt = -i[Htot.Pt]^liVWi]pt^l2V[a2]pt^-fbV[b]pt, (1) 

with V[a]p = apa^ — ^{a^ap + pa^a), which can be 
simplified by transforming to a rotating frame (RFT) 



p't = Uo{typtUo{t) with Uo{t) 



iujt{alai-\-ala2-\-b^ b) 



The RFT leaves the Lindblad terms invariant and gives 

In the small detuning limit, |cjn— ^| "C 2a;, we can neglect 
terms rotating with frequency 2a;, obtaining 



1 TT' 



iao{b - b^) + H,{Aun) + g{Jb^ + J^b), (3) 



where ao = Ad/2h is the Rabi frequency of the driving 
field, J = g~^{g\(7i -\- g20'2) the collective atom-cavity 
interaction, and Hz{Aujn) — Acjicrjcri + Aa;2cr2^2 with 
Acj^ = cjn — cj incorporates the detuning of the atoms 
from the cavity field. The RFT and J above differ slightly 
from the usual interaction picture transformation and 
definition of J to allow for off-resonant and asymmetric 
interactions. We further eliminate zao(6— 6''") by the stan- 
dard displacement transformation p'l = Ub{f^) ptUbif^Y 
with UbiP) = e^^^~^ \ UbiP) commutes with the atomic 
operators an and the standard commutation relations 
[b - 6t,6] = [6,6t] = 1 give Ub{P)aUb{py = a + /3 for 
a = b or b^ (e.g., using Hadamard's lemma), and thus 
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with J^ = J + Jt. Setting P 

term and gives the simphfied master equation 



2 ^ cancels driving field 



p'l = C{p';) - iA76[J&^ + J^b, p'l] + it,V[h\pl (4) 

where the Lindblad operator for the atomic subsystem is 

C{p'l) = -i[H,{Aujn) + aJ,, p'l] + E„ lnV[<Jn]p'l (5) 



and A = g/^bi o^ = Pg = —2aoX. If the cavity is heavily 
damped following the standard procedure for adiabatic 
elimination of the cavity mode (e.g. [3 [10]) leads to the 
equation of motion for the atomic subsystem 



Pt 



Cipt)+jV[J]pt + (V[.P]-V[J]){l\p,\l) (6) 



with 7 = 4A^75 = Ag'^/jb- If the population of the first 
excited cavity mode |1), and thus (l|p^|l), remains small 
and the atomic spontaneous emission rates 7n <C 7 then 
the contributions of {V[J^] -V[J]){l\pt\l) and -fn^[(^n]pt 
will be negligible, yielding the Dicke-type master equa- 
tion for the atomic subsystem with a detuning term 



Pt 



-i[H,{Aoj„) + aJ^,pt] + jV[J]pt. (7) 



If the detuning of the atoms from the cavity frequency 
uj and their interaction with the cavity field are symmet- 
ric, i.e., Auji = Auj = —Auj2 and gi = g2 = g, then 
the detuning term Hz{AuJn) = AcjJ^, and it is easy to 
verify that |^ss) = ^"^(Aa;|n) + a\n) - c^llt)) with 
Q = VAcj^ + 2a^ is an eigenstate with eigenvalue of 
the effective Hamiltonian H = AujJz-\-aJx and the Lind- 
blad operator J and thus a steady state of the dynamics 
([t]). For Auj = (no detuning), the populations of the 
and 1 eigenspaces of J^ are conserved, i.e., the system is 
decomposable and there are infinitely many other steady 
states, rendering |'0ss) non-attractive ^, i.e., if the atoms 
are initially in a state other than the maximally entan- 
gled singlet state (MES) IV^ss), they will not evolve into 
the MES but a low-concurrence mixed state [9], i.e., the 
MES is protected against the decay but cannot be cre- 
ated starting from another state this way. 

To be able to create a MES, or a state arbitrarily close 
to it, we must ensure that the target state is the only 
steady state of the system. This can be accomplished 
by noting that even a small detuning breaks the sym- 
metry and renders l^/^ss) the unique globally attractive 
steady state [5j. Moreover, as the steady state is pure, 
any stochastic quantum trajectory of the system must 
converge to this state. Setting hi = Auj/a^ the concur- 
rence of the steady sate is Cgg = C(^ss) = 2a'^/Q'^ = 
(1 + /t:^/2)~^, showing that attains 1 only for /^ = 0, 
but can be made arbitrarily close to 1 for Auj 7^ by 
increasing a. Fig. [2Va) shows the dependence of Cgs on 
the detuning parameter k,; 99% concurrence is achievable 
with K ^ 0.14. The steady state concurrence Cgg is in- 
dependent of the decay rate 7 but 7 determines the rate 
of convergence to the steady state. Fig. [2|b) shows the 
average time required for the concurrence of the system 
state to reach 99% of Css- This time is virtually indepen- 
dent of the initial state, i.e., the variation for different 
initial states is minimal, but depends strongly on /^, de- 
creasing substantially with increasing k. The graphs also 
suggest that there is an optimum value of 7/a between 
1.5 and 2.5, depending on /^, with j = 2a being close to 
optimal for a wide range of k,. This suggests that the best 




FIG. 2: (Color online): Steady state concurrence C as func- 
tion of K (top) and log^Q of time needed to reach 99% of the 
steady state concurrence (bottom) as function of 7/a and k. 
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FIG. 3: (Color online): Concurrence C{t) (bottom) for vari- 
ous (simple) detuning profiles /\uj{t) (top) shows that expo- 
nential detuning achieves the highest concurrences fastest. 



strategy to achieve high concurrence in the shortest pos- 
sible time is to choose ^/a near the optimum value and 
the largest possible n that still allows us to reach the de- 
sired concurrence, or better yet, start with a large Kq and 
gradually decrease it to 0, e.g., linearly or exponentially. 
To assess whether we can achieve high entanglement 
for the setup in FIG. [l] with this scheme, we turn to sim- 
ulations. For ^^Rb we choose uo = 6.83 GHz, and for 
atoms a few microns apart, detunings up to 1 MHz are 
attainable via a magnetic field gradient. Recent exper- 
imental results [11] suggest that realistic values for the 
coupling constant g and effective Rabi frequency for the 
driving field ao in a lossy microwave cavity are on the or- 
der of a few hundred kHz, and the cavity damping rate 75 
can easily be made about one order of magnitude larger. 
We choose a^ = g = 200 kHz and 75 = 10^, which gives 
A = 0.1, 7 = 80 kHz and a = 40 kHz, and 7/a = 2, 
which is in the optimal range to minimize the conver- 
gence time tgg according to Fig. |2Vb). Our previous re- 
sults show that to achieve a steady state concurrence of 
99% requires n < 0.14 or Acj < 5.6 kHz. Fig. [3] shows 



that with a constant detuning Acj = 5.6 kHz we can 
reach 99% concurrence from an arbitrary initial state in 
less than 20 ms. The figure also shows that this result 
can be significantly improved if we start with a larger 
detuning and gradually decrease it to 0. E.g., starting 
with Acj = 2.5a = 100 kHz, which is still well below the 
limit of 1 MHz, we can obtain over 99% concurrence in 
less than 5 ms if the detuning is decreased exponentially. 

Next we verify the validity of the approximations made 
in deriving ([7|), i.e., adiabatic elimination and negligibil- 
ity of the terms 7^^^]/)^ and {V[J^] - 7:>[J])(l|p^|l). 
Comparison of the exact solution pt given by Q and the 
solution p^ of the approximate model (ItI suggests that 
the latter is a very good approximation in the parame- 
ter regime above. For 71 =72 =0, the reduced density 
operator for the atomic subsystem obtained by solving 
Q and tracing over the cavity modes, p^ = Tr5[p^], and 
the approximate solution p^ of (JTl are in excellent agree- 
ment with the norm error e = max^ ||p^ — p^|| < 0.0201 
for Ao; = 5.6 kHz (const) and < 0.0172 for Aa;(t) = 
100e~^"^^, and we have verified that the concurrences pre- 
dicted by the exact and approximate models are in very 
good agreement for all the simulations below. 

For the systems considered here 7n <C 7 should be eas- 
ily achievable under current experimental conditions, as 
the lifetimes of hyperfine states of ultracold atoms such as 
^'^Rb are on the order of seconds, suggesting 7n < 1 Hz, 
while our effective cavity decay rate is 7 = 80 kHz, but 
the effect of non-zero atomic decay rates on the attain- 
able concurrence deserves investigation. Fig. [4] shows the 
concurrence of the system state as function of time and 
the atomic decay rates 71 = 72 for constant and exponen- 
tially decreasing detuning. For 7^ = 0.001 kHz, the ef- 
fect of the atomic decay is virtually imperceptible on the 
timescales considered but for 7 = 0.1 kHz, atomic decay 
decreases the attainable concurrence substantially. The 
exponentially decaying Acj achieves a higher peak con- 
currence in the presence of atomic decay, but as the de- 
tuning approaches 0, atomic decay takes over and drives 
the concurrence down again, while for constant detun- 
ing the concurrence reaches a lower steady state value 
but remains near constant on the timescales considered. 
This might be explained in terms of inhibition of spon- 
taneous emission by detuning p!2l resulting in reduced 
effective decay rates. This suggests that in the pres- 
ence of non-negligible atomic decay the best strategy is 
to start with a large initial detuning to maximize the 
rate of convergence, but to decrease it not to zero but 
a non-zero asymptotic value Ao;/ (offset), chosen such 
that 2/(2 + (^^)^) is close to the peak value for the 
concurrence. Indeed, Fig.jsjleft) shows that this strategy 
appears highly effective in maximizing the steady-state 
entanglement and the rate at which it is reached. 

Another possible source of error is asymmetry in the 
couplings of the atoms to the cavity. If the atoms are 
identical and interact in phase with the cavity field, i.e.. 




FIG. 4: (Color online): Concurrence C{t) as a function of time 
and atomic decay rates 71 = 72 for Acj = 5.6 kHz constant 
(left) and exponential detuning Acj(t) = 100e~°-^* (right). 
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FIG. 5: (Color online): Concurrence C{t) as a function of 
time and the atomic decay rates 71 = 72 (left) and time and 
asymmetry in the atom-cavity couplings rig (right) for expo- 
nential detuning with constant offset. 



if they are well localized in the cavity, gi = ^2 should be 
a good approximation even in the presence a small mag- 
netic field gradient S as ^ is proportional to the mag- 
netic dipole moment, which is not affected by a small 
magnetic field gradient. If there is an asymmetry, say 
^2 = (1 + %)^i, the system no longer has a pure steady 
state but a unique mixed state in the interior, which has 
lower concurrence than the MES, as illustrated in Fig.|6] 
Again, as for atomic decay, the exponential detuning pro- 
file allows us to reach higher concurrences fast but the 
concurrence decreases as the detuning vanishes. Fig. |5] 
(right) shows that the results can again be substantially 
improved by adding an offset to the exponential detuning 
to stabilize the concurrence near its peak value. 

Finally, we consider the effect of a frequency-offset 
that might result if the atoms are not positioned sym- 
metric with respect to the magnetic field gradient, lead- 








FIG. 6: (Color online): Concurrence C{t) as a function of 
time and asymmetry in the atom-cavity couplings rig for Acj = 
5.6 kHz (left) and Acj(t) = 100e-°-^* (right). 



ing to asymmetric detuning and an error term H^ = 
r]^Auj{alai + 0-2 0-2). Surprisingly high entanglement is 
attainable even for large rj^j . A numerical fit to the simu- 
lation data shows that the concurrence after 20 ms is well 
described by a quadradic model C2o{t]uj) = <^2^^ + <^i^u; + 
ao with (a2,ai,ao) = (-0.062,-0.0024,0.99) for Au = 
5.6 kHz and (-0.15, -0.039, 1.00) for Auj{t) = lOOe-^-^^ 
For r]^ > 0.2 the constant detuning is preferable as the 
large initial amplitude of the exponential detuning mag- 
nifies the error, but in both cases substantial entangle- 
ment, 92.66% and 81.09%, is still attainable even if the 
offset is on the order of the detuning, r]^ = 1. 

We have shown that effectively maximal entanglement 
of non-interacting atoms can be realized by putting the 
atoms in a lossy cavity driven by a classical field and ap- 
plying a magnetic field gradient. While there are many 
schemes for entangling qubits, this scheme exploits a 
novel concept of environment-induced entanglement to 
create high-grade entanglement and offers exceptional 
simplicity^ eliminating the need for direct coupling be- 
tween the atoms, control lasers to manipulate internal 
atomic states, and complicated measurement setups and 
problems due to limited detection efficiency. Global at- 
tractivity of the target state ensures that the atoms con- 
verge to this state for any initial state, eliminating the 
need for initial state preparation. Error analyis shows 
that high-grade entanglement can be still achieved with 
simple detuning profiles in the presence of atomic decay, 
asymmetric atom-cavity coupings or frequency offsets. 
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